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Doppler and ranging measurements between spacecraft can be obtained only 
when the ratio of the total received signal power to noise power density (Pt/No,) 
at the receiving spacecraft is sufficiently iarge that reliable signal detection can 
be achieved within a reasonable time period. In this article, the requirement on 
Pt/No for reliable carrier signal detection is calculated as a function of various 
system parameters, including characteristics of the spacecraft computing hardware 
and a priori uncertainty in spacecraft-spacecraft relative velocity and acceleration. 
Also calculated is the Pt/No requirement for reliable detection of a ranging signal , 
consisting of a carrier with pseudonoise (PN) phase modulation . Once the Pt/No 
requirement is determined, then for a given set of assumed spacecraft telecommuni- 
cation characteristics (transmitted signal power, antenna gains, and receiver noise 
temperatures) it is possible to calculate the maximum range at which a carrier 
signal or ranging signal may be acquired. For example, if a Mars lander and a 
spacecraft approaching Mars are each equipped with J-m-diam antennas, the trans- 
mitted power is 5 W, and the receiver noise temperatures are 350 K, then S-band 
carrier signal acquisition can be achieved at ranges exceeding 10 million km. An 
error covariance analysis illustrates the utility of in situ Doppler and ranging mea- 
surements for Mars approach navigation. Covariance analysis results indicate that 
navigation accuracies of a few km can be achieved with either data type. The analy- 
sis also illustrates dependency of the achievable accuracy on the approach trajectory 
velocity. 
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I. Introduction 

In recent years, there has been increased interest in pur- 
suing an intensive program of lunar and Mars exploration. 
Conducting an expansive program of lunar and Mars ex- 
ploration with multiple spacecraft will present new tech- 
nical and operational challenges. For example, future mis- 
sions to Mars may rely on aerobraking to achieve orbit in- 
sertion. The aerobraking technique reduces a spacecraft’s 
propellent requirement, but imposes more stringent navi- 
gation constraints than have been needed for previous in- 
terplanetary missions. Simultaneous deployment of a large 
number of spacecraft at the Moon and/or Mars could place 
unprecedented demands on already taxed ground-based 
tracking resources. The need for greater navigational ac- 
curacy and the need to relieve the burden on ground-based 
antennas could both be addressed by making use of in situ 
radio-metric tracking measurements. Doppler and ranging 
measurements between spacecraft provide navigation in- 
formation complementary to Earth-based techniques and 
could potentially decrease reliance on ground-based track- 
ing. 

This article investigates requirements for successful ac- 
quisition of in situ radio-metric data types and discusses 
potential measurement precision. An error covariance 
analysis has also been performed to investigate the nav- 
igation accuracy that can be achieved using Doppler and 
ranging measurements between a spacecraft approaching 
Mars and a Mars-orbiting satellite used as a navigation 
aid. The analysis was designed to build upon the results 
of similar studies conducted previously [1-4]. The error 
covariance analysis points out some of the most important 
trajectory and tracking system parameters that influence 
the navigation performance that can be obtained. While 
the mission scenario studied herein involves Mars, in situ 
radio tracking is also applicable to lunar spacecraft navi- 
gation, or to any other mission in which a radio beacon or 
another spacecraft is available near the target body and 
can be used as a navigation aid. 

II. Carrier Signal Acquisition 

To perform in situ Doppler measurements, one space- 
craft tracks the phase of a carrier signal transmitted from 
a second spacecraft. If the measurement is one-way, the 
signal originates at the second spacecraft; in a two-way 
measurement, the signal is uplinked by the first spacecraft, 
coherently transponded by the second spacecraft, and then 
phase tracked by the first spacecraft. The Doppler shift 
on the received signal provides a direct measure of the rel- 
ative line-of-sight, velocity of the two spacecraft. While 
one-way measurements are more easily implemented, they 


are inferior to two-way measurements because an unknown 
frequency offset between the two spacecraft oscillators can 
cause a significant error in the inferred relative velocity of 
the spacecraft. 

A. Detection Strategy 

A simple signal detection strategy is depicted in Fig. 1. 
The incoming signal is mixed with a locally generated 
model; the resulting baseband signal is integrated for T 
sec; and N consecutive T - sec integrations are Fourier 
transformed. The resulting Fourier coefficients reflect sig- 
nal or noise amplitude in N bins of width 1/(NT) Hz, cen- 
tered at the model frequency (Fig. 2). In general, the a pri- 
ori uncertainty of the signal frequency, Av, will be larger 
than the frequency range, 1/T Hz, that can be searched 
with a single Fourier transform, so a number of different 
model frequencies will have to be tried before the signal is 
detected. 

To determine whether the signal is present in a given 
frequency span (the range of frequencies probed by a single 
Fourier transform), an amplitude detection threshold is 
selected. It is desirable to have the threshold be sufficiently 
low so that the signal will almost always be detected when 
it is present, and sufficiently high so that noise will not 
often be mistaken for signal. If the amplitude in some 
Fourier resolution bin exceeds the threshold, then the same 
frequency model is used to process the next NT sec of data. 
If a signal is detected in the same bin in two consecutive 
Fourier transforms, then it is assumed that the true signal 
has been found, and phase-lock-loop tracking commences. 

Note that in any given Fourier transform, it is possible 
that the noise amplitude in at least one Fourier resolution 
bin will exceed the amplitude detection threshold. If one 
requires 99-percent confidence of detecting the true sig- 
nal twice consecutively when it is present, and 1- percent 
probability of mistaking noise for signal, one obtains a re- 
lation between N and ( P t /No)T that is approximated by: 

( P t /No)T = 29N~ 0 92 . The Appendix discusses how this 
equation is obtained. 

If the signal is equally likely to lie anywhere within the 
search space, A*/, then the mean time to signal acquisi- 
tion is T acq & At/NT 2 /2 [5]. The acquisition time is 
proportional to NT because that is the time required to 
accumulate data for each Fourier transform, and is propor- 
tional to T Az/ because the frequency search space contains 
TAi / spans of width (l/T) Hz. 

Taken together, the equations ( P t /No)T — 29N~ 092 
and T AC q » AvNT 2 /2 imply that for a given value of 
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either N,T, or NT, there is a functional relationship be- 
tween Tacq/Av and the minimum P t /No required for re- 
liable signal acquisition (99-percent confidence of signal 
detection and 1-percent probability of false alarm). The 
contours in Fig. 3 illustrate these relationships. As an 
example, if N = 1000 and T = 0.01 sec, it is required 
that (Pt/No) > 7 dB-Hz. The corresponding value of 
TacqI Ai/ is about 0.05 sec/Hz; e.g., for a frequency un- 
certainty of 1 kHz, the expected search duration is 50 sec. 

It is important to realize that the quantities N and T 
are not free parameters. Constraints due to both space- 
craft dynamics and hardware (and perhaps software) im- 
pose restrictions on the allowable values of N, T, and their 
product, NT. 

1. Constraints on the Selection of N and T. Hard- 
ware on board the spacecraft dictates a maximum value 
of N and a minimum value for T. At present, Fourier 
transform chips that can operate on hundreds of points 
are commonplace. If N is constrained by hardware to be 
at most 1024, then the data dump interval, 7\ is greater 
than 10“ 4 sec over most of the range of Pt/No in Fig. 3. A 
dump rate, 1/T, of 10,000 Hz is well within the capabilities 
of existing ground-based hardware and may be feasible for 
space-qualified systems by the late 1990’s. 

2. Constraints on Data Collection Interval, NT 

Constraints on the product NT result from unmodelled 
beacon accelerations and from oscillator drift, both of 
which cause the frequency, /, of the received signal to 
change in time. Signal detection can be severely impeded 
if the signal moves through many Fourier resolution bins 
over the time, NT, during which the data for a single 
Fourier transform are acquired. If it is required that the 
change in frequency due to an unmodelled acceleration, 8a , 
over time NT be smaller than l/(2iVT) Hz (half the size 
of one Fourier resolution bin), then the following relation 
results: 


NT ^'{w« m 

A change in frequency due to oscillator instability has 
an effect similar to an unmodelled acceleration. Let 
<r yi (NT) be the Allan standard deviation of spacecraft Fs 
oscillator. Then imposing a constraint analogous to the 
unmodelled acceleration constraint above yields 

NT < 1 = (2) 

2/^, (A rT) + al(NT) 


Contours corresponding to NT = 10 n for n — —2 to 2 
are plotted in Fig. 3. 

3. Range of Frequencies to be Searched, A v. The 

range of frequencies to be searched, A v, is a function of 
the uncertainty in spacecraft -spacecraft relative velocity. 
Also, for one-way transmission (and for the uplink of a two- 
way measurement), Nv depends on an unknown frequency 
offset in the transmitting spacecraft’s reference oscillator. 

If the spacecraft-spacecraft velocity uncertainty is 8v, 
the corresponding uncertainty in carrier frequency is Nv - 
f(8v/c). The day-to-day drift in carrier frequency due 
to oscillator instability depends on the quality of the 
spacecraft oscillator. If the oscillator is not temperature- 
stabilized, a 2.3-GHz (S-band) carrier frequency may vary 
by thousands of Hz. A highly stable oscillator has day-to- 
day variations much smaller than I Hz at S-band. 

4. Example: Calculation of Pf/N 0 Requirement 
for Carrier Acquisition. The above information may 
be combined in order to calculate the Pt/No requirement 
for reliable carrier acquisition. 

As an example, suppose one spacecraft is fixed on the 
surface of Mars, while a second spacecraft is on approach 
to Mars. The spacecraft are equipped to transmit and 
receive a 2.3-GHz signal. It is desired that a carrier sig- 
nal be acquired within 5 min. Neither spacecraft has a 
priori knowledge of the spacecraft -spacecraft relative ve- 
locity and acceleration. 

Since the spacecraft have no a priori knowledge of their 
relative velocity, the range of frequencies to be searched, 
A v, is the full spacecraft-spacecraft relative velocity. A 
typical Mars approach velocity is 4 km/sec, corresponding 
to a frequency search range of approximately 31 kHz at 
S-band. Since 5- min acquisition is desired, it is required 
that T acq! ~ 0.01 sec/Hz. 

The relative acceleration of the lander and approach 
spacecraft is principally due to Mars rotation, which im- 
parts up to 0.017 m/sec 2 of acceleration to the lander. This 
unmodelled acceleration imposes the constraint NT < 
2 sec. Equation (2) indicates that the upper bound on 
NT arising from oscillator stability is also 2 sec if the 2- 
sec Allan standard deviation of the spacecraft oscillators 
is smaller than 8 x 10 -11 sec/sec. 

The parameter constraints are summarized in Table 1. 
The constraints TacqI Ai' ~ 0.01 sec/Hz and NT < 2 sec 
are represented by the unshaded region in Fig. 4. The min- 
imum Pt/No satisfying both constraints is approximately 
13 dB-Hz. 
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Once the Pt/No requirement is known, then if the 
spacecraft telecommunication characteristics are specified, 
one can determine the maximum range at which the carrier 
can be acquired. Figure 5 illustrates P t / No as a function 
of range between the Mars lander and approach space- 
craft, for several different assumptions about the space- 
craft telecommunications system. The tapering of the 
curves at smaller ranges is due to the increasing contri- 
bution of Mars to the system noise temperature. 

III. Doppler Measurement Precision 

Once a carrier signal has been detected, tracking the 
signal’s phase is straightforward. Phase tracking precision 
depends on a variety of factors, including Pt/No and os- 
cillator stability. However, if the carrier can be tracked at 
all, then the point-to-point phase precision is much smaller 
than one RF cycle (e.g., 13 cm at S-band). Phase tracking 
precision does not significantly affect the accuracy with 
which the Mars approach spacecraft’s state can be deter- 
mined by fitting to an arc of beacon-spacecraft Doppler 
data. Determination of the approach spacecraft’s state 
should be limited by other effects, such as unmodelled 
spacecraft accelerations. 

IV. Ranging Signal Acquisition 

A. Signal Structure 

The ranging measurement is performed as follows: A 
spacecraft transmits a signal to a second spacecraft that 
locks onto the signal and coherently transponds it back to 
the first spacecraft. Two-way ranging offers the tremen- 
dous advantage of eliminating the problem of unknown 
spacecraft clock offsets that can significantly degrade the 
accuracy of a one-way measurement. However, the accu- 
racy of the two-way measurement is still affected by the 
first spacecraft’s clock instability over the round-trip light 
time. For example, if the spacecraft oscillator has an Allan 
variance of 10“ 11 sec/sec over time scales of minutes and 
the two-way light time is 100 sec (30 million km), then 
clock instability introduces a 30-cm error. 

The ranging signal is assumed to be of the following 
form: 

s(|) = A cos (u; 0 f -f ftP N (p, f c ,t -* n ) + 0 o ) 

where 


u>q = the carrier frequency 
(3 = the modulation index, 0 < ft < n/2 

PN = a pseudonoise sequence with period p and chip 
rate f c 

p — the period of the PN sequence 

f c = the chip rate of the PN sequence, i.e., 1 / f c is the 
chip duration, r c 

t n = the epoch of the PN sequence 

0o = the phase of the carrier signal at t = 0 

Pseudonoise (PN) sequences [6] are a special class of 
sequences of l’s and — T’s that can be generated by linear 
feedback shift registers. Individual elements of such a se- 
quence will be referred to as “symbols.” PN sequences are 
periodic; a given sequence consists of repetitions of a sym- 
bol pattern of length p symbols. For every integer n, there 
exists at least one PN sequence of period 2 n — 1. If r(j ) is 
a PN sequence with period p, then the autocorrelation of 
r(j),Rp(i) = (1 fp) ff r(j)r( j -f 1 )dj is equal to 1 for i = 
any integer multiple of p and is equal to -1/p for all other 
i (Fig. 6). The chip rate, / c , of a PN sequence is the rate 
at which symbol transitions occur. 

During a two-way ranging measurement, the spacecraft 
originating the PN signal simultaneously receives a PN 
code that had been transmitted one round-trip light time 
earlier. The incoming PN code is offset from the current 
state of the uplink code by an amount equal to the round- 
trip light time divided by the chip duration, r c . This code 
offset can be determined by cross-correlating the received 
signal with the locally generated sequence. To ensure that 
such a two-way ranging measurement is unambiguous, it is 
required that the period of the PN sequence, pr c , be larger 
than the a priori uncertainty in the spacecraft-spacecraft 
range. The relationship between the chip rate, PN period, 
and range ambiguity is plotted in Fig. 7. The ranging 
precision is proportional to the product of the chip dura- 
tion and the precision with which the PN sequences can 
be aligned. If chip alignment precision of 1 percent can be 
achieved in the tracking loop, then a chip rate of 1 MHz 
corresponds to a ranging measurement precision of 3 m. 

B. Ranging Signal Detection 

The incoming signal is 


.4 = the signal amplitude s (t) = A cos (w 0 < + 0PN(p, f e , t - t„) + <f> 0 ) 
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and the locally generated model for the carrier is 
s(t) = exp [i(&ot + 0PN(p, fcJ- In ) )] 


If the carrier is fully suppressed (0 — tt/ 2), then mul- 
tiplication of the signal and model, followed by low-pass 
filtering, yields the following: 


C(0 = s(t)(s(t)) = j exp 


i ((lD 0 — LOo)t — (po) 


Even if the model frequency is perfect, a fully sup- 
pressed carrier cannot be detected unless the model PN 
epoch is also close to truth (within a fraction of a chip). If 
the model PN epoch is stepped in single-chip increments, 
then the best PN code alignment achieved may be off by 
as much as 1/2 chip, resulting in a factor of two ampli- 
tude loss (see Fig. 6). In order to decrease the potential 
loss of amplitude due to fractional PN-code misalignment, 
one may conduct the search using smaller step sizes, at 
the cost of requiring additional models to cover the entire 
search space. If the step size is 1 jk c chip, then up to k c x p 
offsets may need to be tested before the signal is detected. 


x 


PN(pJ c ,t 


tn)PN(pJct-in) 


To detect the signal, it is necessary that the model pa- 
rameters, w 0 , /c and i n , be adequately close to the actual 
signal parameters, wo,/ c , and t n . However, it is possi- 
ble to design the spacecraft so that the PN clock and the 
transmitted carrier frequency are both referenced to the 
same frequency standard. In this case, a search must be 
conducted only for w 0 and t n \ the code frequency, / c , is a 
deterministic function of the carrier frequency, wq. 

The product of actual and model signals is integrated 
over T - sec intervals, and a set of N consecutive inte- 
grations is Fourier transformed. It is important that 
the time period NT be long enough that the product 
PN(p,f c ,t — in) [PN{p,fc,t - /„)] approximates the PN 
autocorrelation function. The signal can be detected \ia 
the Fourier transform only when (l) the model frequency 
lies within 1/27' IIz of t he actual frequency and (2) the 
difference between the model PN epoch and the true PN 
epoch is a fraction of a chip. 

Since a given Fourier transform probes a region of fre- 
quency space of width 1// Hz, centered at ihe model he- 
quency. a total of VAiz Fourier transforms are required to 
search a frequency range of A/z Hz. Note, however, that a 
7-sec average of {signal x model} suffers an amplitude 
loss equal to 1 * * sine (7 1 (w () — wo/2); il the signal and model 
differ by 1/27’ Hz. the corresponding loss in amplitude is 
about 40 percent. It may be advantageous, therefore, to 
space the frequency models by an amount \/kjT{kj > 1) 
to reduce the distance between the signal and the best 
model frequency. Reducing the sparing between frequency 
models increases the number of models needed to cover the 
search space. 

1 S. A. Stephens. “An Analysis of FI T Tone Acquisition. JIT In- 

teroffice Memorandum TFT 1-92- I t (internal document). Jet Pro- 

pulsion Laboratory, Pasadena, California. May 11, 1992. 


The signal search space is illustrated in f ig. 8 for the 
case kj = 1. The search space consists of Nj x N s cells, 
where N f is the number of model frequencies and N s is the 
number of trial symbol offsets. Each Fourier transform of 
{signal x model } probes a range of (1 /T) Hz, centered 
at the model frequency. One method of performing a sys- 
tematic search is to step through all possible model symbol 
offsets for a given model frequency. If the signal is not de- 
tected, the model frequency is then shifted by 1/7 Hz. 
Thus, at each model frequency, N s Appoint, Fourier trans- 
forms are performed. Assuming the Fourier transforms can 
be computed in real time, N s x(NT) sec are spent search- 
ing at each model frequency. Unfortunately, it is possible 
that this search algorithm would miss t he signal entirely if 
the incoming signal's frequency is changing with time (due 
to u n modelled spacecraft- spacecraft relative acceleration 
and oscillator drift). 

The problem of a signal with a moving frequency may 
be addressed by taking the following two steps: (1) set k f 
equal to 2, i.e. , space the model frequencies by 1/2T Hz; 
each Fourier transform of {signal x model) will still probe 
a region in frequency space of width 1 /T Hz, but the region 
probed by adjacent frequency models will now overlap by 
ot) percent; and (2) require that, the change in signal fre- 
quency over a time \ s x ( A / ) be smaller than (1/27 ) Hz. 

( Jive n an upper bound on the spacecraft spacecraft, rel- 
ative acceleration and oscillator drift, this requirement- 
places a constraint on the allowable values of the parame- 
ters A7, N , and T (the impact of this and other constraints 
is discussed later). The number of Fourier transforms re- 
quired to search the entire two-dimensional space with this 
st rategy is 27 A v AN . 

The signal detection procedure is the same as the car- 
rier detection method discussed earlier. Model signals are 
sequentially tested by mixing with tin' incoming signal and 
then Fourier transforming A 7 points acquired over A / sec. 
If a Fourier frequency bin has amplitude exceeding the de- 
tection threshold, then the same model is tested again. If 
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a second detection occurs, it is assumed that the signal 
has been found and phase- and delay-lock-loop tracking 
is commenced. The mean time to signal acquisition is 
T acq & NT 2 Ai/N 9 . The mean time to acquisition could 
be significantly reduced if the spacecraft hardware is ca- 
pable of processing multiple model signals in parallel. 


that the range does not change by more than 1/2 chip 
(r c /2 sec) over the NT sec during which data are collected 
for a Fourier transform. 


NT < 


2 Sv 


(4) 


If the signal power is known, the amplitude threshold 
can be selected so that the probability of signal detection 
when a signal is present, Pdetect , is high (e.g., 99 per- 
cent), while the probability of a false alarm, Pf a i se , is 
low (e.g., 1 percent). For a given probability of detec- 
tion and probability of false alarm, there is a functional 
relationship between the quantities (P t /N$)T and N. For 
Pdetect = 0.99, Pjaise — 0.01, and k c — 2 (symbol search 
proceeding in half-chip steps), the relationship between N 
and (P t /N 0 )T is well-approximated by the following equa- 
tion (see Appendix): 

Aj' = 43 ,v -0 91 
/Vu 

The contours in Fig. 9 are derived by combining the 
equations T AC q = NT 2 AvN s and (P t /N 0 )T = 43 A" 091 . 
As in the case of pure carrier detection, there are a number 
of constraints that restrict the allowable values of A r ,T, 
and NT. 

1. Constraint on NT Arising From Unmod- 

elled Spacecraft Acceleration. The requirement that 
the change in signal frequency over a time N s x (NT) be 
smaller than (1/2T) Hz is expressed as follows: 

■m < m. sr- < _ f_ 

dt N S NT 2f 0 Ng6a 

Here f 0 is the carrier frequency, c is the speed of light., and 
ba is the unmodelled spacecraft acceleration. As in the 
case of simple carrier detection, it is also required that the 
time A T be small enough so that mimodelled spacecraft 
accelerations do not smear the signal over multiple Fourier 
frequency bins. Together, these two constraints are given 
by 

2. Constraint on NT Arising From Unmodelled 
Spacecraft Velocity. It is required that the unmodclled 
spacecraft -spacecraft relative velocity be small enough so 


where 6v is the unmodelled spacecraft-spacecraft relative 
velocity. 

3. Oscillator Stability Requirements. The con- 
straints on oscillator stability are analogous to the con- 
straints on NT arising from unmodelled spacecraft accel- 
eration and velocity. 

In analogy with the acceleration constraint, it is re- 
quired that the signal change by less than (1/2T) Hz over 
N s x (NT) sec, and also that the carrier frequency not 
move through many Fourier frequency bins over the NT 
sec during which data are acquired for the Fourier trans- 
form (A frequency < 1 /(2NT) over NT sec). Letting 
o'yfr) denote the r-sec Allan standard deviation of space- 
craft f s oscillator, 

VT) < AA 

(5) 

\[T- (NT) + <rl (NT) < 

V Sl ' ’ 'J- y ’ 2 f 0 NT 

In analogy with the velocity constraint, it is required 
that oscillator drift not introduce a delay of more than 
r c f 2 sec over the NT sec during which data are acquired 
for a Fourier transform. This constraint is looser than that 
imposed by the above equations for chip rates smaller than 
the carrier frequency (the usual case). 

4. Example: Calculation of P t /N 0 Required for 
Ranging Signal Acquisition. Suppose that ranging 
measurements are to be performed between a spacecraft 
approaching Mars and a beacon spacecraft on the planets 
surface. The carrier frequency is S-band, and the chip- 
ping rate is 1 MHz. The length of the PN code is 2047 
symbols, corresponding to a two-way range ambiguity of 
600 km. The ambiguity can be resolved via spacecraft po- 
sition determinations resulting from routine Earth-based 
Doppler tracking of the two spacecraft. 

The spacecraft are provided with a crude estimate of 
their relative velocity: 0.5 km/sec, equivalent to 4000 Hz 
at S-band. The spacecraft-spacecraft relative accelera- 
tion is not modelled. Until the approach spacecraft is 
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within the sphere of influence of Mars (£1 million km), 
the total relative acceleration of the spacecraft and bea- 
con is 0.017 m/sec 2 due to Mars rotation. Finally, it is 
assumed that the spacecraft hardware is capable of per- 
forming 1024-point. Fourier transforms. 

The parameter constraints corresponding to a signal ac- 
quisition time of 5 min are summarized in Table 2. All four 
constraints are satisfied in the unshaded region of Fig. 10 
bounded by the contours N = 1024 and Tacq/2N s Ai/ — 
10~ 5 sec/Hz. The minimum Pt/No satisfying the con- 
straints is approximately 28 dB-Hz. Equation (5) may 
now be applied to determine the oscillator stability re- 
quirements: 


<r y (409 sec) < 2 x 10~ 6 sec/sec 

cr^O.l sec) < 2 x 10~ 9 sec/sec 

These Allan variances are easily met by non-temperature- 
stabilized crystal oscillators. 

Figure 5 indicates that if both the lander and ap- 
proach spacecraft are equipped with 1-m-diam antennas 
and mediocre receivers ( T rec = 350 A'), a 5-W S-band 
signal is sufficient to enable a 28-dB-Hz signal power-to- 
noise density ratio at spacecraft separations of millions of 
kilometers. 


V. Ranging Measurement Precision 

To within a factor of order unity, the precision with 
which PN codes can be aligned with a delay lock loop 
is l/y/2P t /N 0 T dweU , where Tdweii is the integration time. 
The resulting ranging measurement precision is 


e p (meters) 


sf^Pt/NvTdwell 


The exact expression for the ranging measurement pre- 
cision depends on details of the delay-lock-loop design. In 
practice, due to dispersion effects in the processing signal 
path, the best alignment precision that can be achieved is 
between 0.01 and 0.001 chips. 

The utility of ranging is to eliminate the unknown bias 
inherent in carrier phase (Doppler) measurements. Once 
this bias has been determined, range changes can be pre- 
cisely tracked by measuring the carrier phase. The phase 
of the carrier can be recovered with precision of a fraction 
of an RF wavelength, even when the transmitted signal 
has a modulation index of tt/2 (fully suppressed carrier). 


VI. Error Covariance Analysis for Mars 
Approach Navigation 

In this section, the navigational utility of in situ radio 
tracking is illustrated through some examples drawn from 
a hypothetical future Mars mission, in which a spacecraft 
approaching Mars acquires Doppler and ranging measure- 
ments from a Mars-orbiting communications relay satel- 
lite. Error covariance computations for a sequential esti- 
mation scheme were performed using approximate models 
of the approach spacecraft and relay satellite trajectories, 
and the Doppler and ranging measurement errors. This 
analysis was not designed to be an exhaustive treatment 
of the approach navigation problem, but merely to illus- 
trate some of the more significant aspects of the use of 
in situ Doppler and ranging measurements in this type of 
application. It should be remembered that although the 
particular application discussed below refers to Mars, in 
situ radio tracking is also applicable to missions to the 
Moon and other target bodies where radio beacons or ad- 
ditional spacecraft may be used as navigation aids. 

The error analysis was performed using the model pa- 
rameters summarized in Table 3. The approach spacecraft 
was assumed to move on a hyperbolic path with respect 
to Mars; except for the final hours prior to arrival at peri- 
apsis, this assumption leads to a trajectory that is essen- 
tially rectilinear and of constant velocity. To illustrate the 
effect of different approach velocities on navigation per- 
formance, two different values of asymptotic approach ve- 
locity were considered, representative of the minimum and 
maximum possible values for low-energy ballistic Earth-to- 
Mars transfer trajectories. The radius of closest approach 
chosen for this study (3417 km) yields a periapsis altitude 
of about 20 km, which is representative of t rajectories lead- 
ing to a direct entry and landing, or the initiation of an aer- 
obraking maneuver for orbit insertion. The Mars-orbiting 
relay satellite was assumed to be in a circular orbit with 
a 12-hr period. The spacecraft acceleration-process noise 
level used in this analysis is intended to account for the 
small, nongravitation al forces (e.g., solar radiation pres- 
sure mismodelling, gas leaks from valves, and pressurized 
tanks) that act on an interplanetary spacecraft. 

The navigational utility of both one-way integrated 
Doppler measurements and two-way ranging measure- 
ments has been investigated. Integrated Doppler data are 
effectively a measurement of the accumulative change in 
the spacecraft-to-spacecraft range over the contact period. 
As shown by the parameters given in Table 3, integrated 
Doppler provides a highly precise measurement of range, 
but with a large constant bias. The frequency drift uncer- 
tainty assumed for the Doppler measurements is represen- 
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tative of oscillator stability on the order of 10” 12 sec/sec. 
Two-way ranging measurements are much less precise than 
Doppler measurements, but provide an accurate measure 
of the spacecraft-spacecraft separation (Table 3). As de- 
scribed earlier, acquisition of a ranging signal requires a 
more extensive search than simple carrier signal acquisi- 
tion. For fixed assumptions about the spacecraft onboard 
computing resources and the desired mean time to signal 
acquisition, the requirement on Pt/No is higher to support 
ranging acquisition. As a result, the acquisition range for 
ranging data is considerably smaller than for Doppler data. 

The results of the error covariance analysis are sum- 
marized in Figs. 11 and 12. These figures show the un- 
certainties in the approach spacecraft encounter coordi- 
nates as a function of time prior to arrival, expressed 
in an aiming plane (B-plane) coordinate system. 2 Fig- 
ure 11 shows the semimajor axis of the aiming plane dis- 
persion ellipse and the linearized time-of- flight uncertainty 
(position uncertainty in the S unit direction divided by 
r^) for an approach velocity of 3 km/sec, while Fig. 12 
shows the same two quantities for an approach velocity of 
6 km/sec. (Mars Observer’s approach velocity was about 
2.4 km/sec; MESUR Pathfinder’s approach velocity will be 
about 5.5 km/sec.) Note that, in practice, delay-lock-loop 
tracking of a PN-modulated carrier signal enables recov- 
ery of the carrier phase, thus enabling simultaneous acqui- 
sition of ranging and integrated Doppler measurements. 
However, Figs. 11 and 12 show the navigational uncer- 
tainties resulting from the use of either data type inde- 
pendently, in order to illustrate differences in the abilities 
of the two data types to determine various components of 
the approach trajectory. 

Both Figs. 11 and 12 illustrate a clear difference be- 
tween Doppler and ranging performance. Since rang- 
ing data provide a direct measurement of the spacecraft- 
spacecraft separation, the time of flight (which in this case 
is essentially the spacecraft range) can be very accurately 
determined immediately upon acquisition. In contrast, 
Doppler data can only measure this component of the flight 
path indirectly and need to be collected for up to several 
days before sufficient information is obtained. This occurs 
because Doppler data depend upon the relative acceler- 
ation of the approach spacecraft and the relay satellite 
(which is inversely dependent upon their relative range) 


2 The aiming plane, or £?-plane, coordinate system is defined by three 
unit vectors, 5, T, and R. S is parallel to the incoming asymptote 
of the approach hyperbola; T is parallel to the Martian equatorial 
plane; and R completes an orthogonal triad with S and T. The 
aim point for a planetary encounter is defined by the miss vector, 
B, which lies in the T-R plane, and specifies where the point of 
closest approach would be if the target planet had no mass and did 
not deflect the spacecraft’s flight path. 


to determine the time of flight. Additional experimenta- 
tion found that the approach-spacecraft acceleration pro- 
cess noise was the principal factor limiting the integrated 
Doppler’s ability to determine the time of flight. 

Figure 12 shows that when the approach velocity is rel- 
atively high, the smaller acquisition distance associated 
with two-way ranging results in the aiming plane disper- 
sion ellipse of the approach spacecraft remaining relatively 
large until about 15 hr prior to arrival. If a final maneu- 
ver is needed to correct the spacecraft’s aim point, this 
implies that the maneuver may have to be performed very 
near the arrival point, with the ensuing risk that there may 
not be sufficient time remaining to redetermine the trajec- 
tory and correct any maneuver execution errors that might 
have occurred. Thus, if large approach velocities must be 
accommodated, early acquisition of Doppler followed by 
joint acquisition of Doppler and ranging data would be 
desirable. 

VII. Summary 

The Pi/Nq required for acquisition of carrier or rang- 
ing signals between two spacecraft is a function of a num- 
ber of parameters, including a priori uncertainty in the 
spacecraft-spacecraft relative state, and the capabilities 
of the spacecraft processing hardware and software (e.g., 
the maximum number of points that can be Fourier trans- 
formed in real time). Once these parameters are speci- 
fied for a given set of assumptions about the spacecraft 
telecommunications systems, one may calculate the maxi- 
mum range at which a signal can be acquired. As an exam- 
ple, it has been shown that S-band radio-metric measure- 
ments between a Mars lander and Mars approach space- 
craft equipped with 1-m-diam antennas are feasible over 
distances of millions of kilometers. The analysis presented 
here can easily be applied to a wide variety of other mission 
scenarios. 

To illustrate the utility of in situ radio-metric measure- 
ments between spacecraft, a covariance analysis was per- 
formed for the case of Doppler and ranging measurements 
between a Mars orbiter and a spacecraft on Mars approach. 
The results obtained are consistent with previous studies, 
and suggest that approach navigation accuracies of a few 
kilometers at Mars may be obtained with either in situ 
Doppler or ranging measurements. 

In the near-term, in situ ranging and Doppler data 
between spacecraft would likely be relayed to Earth for 
navigation processing. Ultimately, onboard computing re- 
sources could be used to perform the navigation updates, 
enabling a near-real-time in situ navigation capability. 
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Table 1. Parameter constraints for 5-min carrier 
signal acquisition between a Mars lander and a 
spacecraft on Mars approach. 


Constraint 

Reason 

T aCq/^v ~ 1°“ 2 sec/Hz 

5-min acquisition, = 31 kHz 

NT < 2 sec 

Eq. (1) (unmodelled acceleration) 


Table 2. Parameter constraints for 5-min ranging signal acquisi- 
tion between a Mars lander and a spacecraft on Mars approach. 

Constraint 

Reason 

T ACq! 2N s Au £ 10 -5 sec/Hz 

5-min acquisition, A^ = 4000 Hz, 
and N s = 4094 

N < 1024 

Capability of Fourier transform chip on 
spacecraft 

NT < 0.03 y/N sec 

Eq. (3) (unmodelled acceleration) 

NT < 0.3 sec 

Eq. (4) (unmodelled velocity) 
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Table 3. Mars approach navigation error analysis parameters. 


Trajectory parameter 

Value 

Asymptotic approach velocity, v oo 

3 or 6 km/sec 

Radius of closest approach 

3417 km 

Inclination of approach trajectory to Martian equator 

+ 15 deg 

Mars relay satellite orbital period 

2 hr 

Mars relay satellite inclination 

0 deg 

Spacecraft acceleration process noise 
(approach spacecraft and relay satellite) 

1.4 X 10~ 5 m/sec/\/hr 

A priori state uncertainty (lcr) 

Value 

Approach spacecraft epoch position 

100 km 

Approach spacecraft epoch velocity 

1 m/sec 

Relay satellite epoch position 

3 km 

Relay satellite epoch velocity 

0.1 m/sec 

Doppler measurement parameter 

Value 

Maximum acquisition range 

2 X 10 6 km 

Additive random measurement noise 

4 mm (lcr) 

Accumulative random measurement noise 

0.58 m 2 /hr 

Frequency drift stability 

0.3 mm/sec (la, a priori) 

Carrier phase bias uncertainty 

100 km (1<7, a priori) 

Ranging measurement parameter 

Value 

Maximum acquisition range 

8 X 10 5 km 

Additive random measurement noise 

5 m (lcr) 

Range bias uncertainty 

15 m (lcr) 
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Fig. 1. A simple carrier signal detection strategy. 
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Fig. 2. The relationship between frequency search range, 

Hz, and frequency span, 1/7 Hz, probed by a single Fourier 
transform. 
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Fig. 3. Relationships between T ACQ lAv, N, 7 , jjnd the minimum 
P t /Nq required for carrier signal acquisition with P fa/se - 0.01 
while P detect- 0.99. 




Fig. 7. Relationship between chip rate, length of PN sequence, 
range ambiguity, and ranging measurement precision (assumed 
tracking error is 1 percent of chip length). 


Fig. 4. Given the constraints NT < 2 sec and T ACQ Mv < 0.01 
sec/Hz, the minimum P t IN 0 needed for reliable carrier detection is 
13 dB-Hz. 
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Fig. 8. Search space for carrier signal with PN phase 
modulation. Each bracket pair represents the range of 
frequencies probed by a Fourier transform of N 7-sec signal x 
model points. 


Fig. 5. P t fN 0 as a function of Mars lander-approach spacecraft 
range. Transmitted power is 5 W. T rec is the receiver noise 
temperature. 
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Fig. 9. Relationships between 7^ C Q/(2W s 4vj, N, T, and the P t 
t N 0 required for ranging signal acquisition with Pf a i se = 0 01 
while P de tect = 0.99. 



Fig. 10. Given the constraints T acq H2N s Av) < 10~ 5 and N < 
1024, the minimum PJNq enabling reliable ranging signal 
detection is 29 db-Hz. 




Fig. 11. Mars approach navigation dispersions ( 3 km/sec): (a) aim-point uncertainty and (b) linearized time-of-f light uncertainty. 




Fig. 12. Mars approach navigation dispersions ( = 6 km/sec): (a) aim-point uncertainty and (b) linearized time-of-flight uncertainty. 
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Appendix 

Calculation of Relations between JVand (P t /N 0 )T 


As described in the article, the carrier detection strat- 
egy involves looking for Fourier resolution bins whose am- 
plitude exceeds a threshold level. The detection threshold 
parameter, £, is selected such that the probability that the 
amplitude of the signal plus noise exceeds 
£ x {expected signal amplitude } is 0.995. The probabil- 
ity of two consecutive signal detections when the signal is 
present, which shall be referred to as Pdetect , is then 0.99. 
Calculation of the detection threshold must take into ac- 
count signal amplitude losses resulting from two effects: 
(1) If the signal lies within the l/T Hz range probed by 
the Fourier transform, it may differ from the model by up 
to l/T Hz. Thus, each T - sec average of the signal x model 
phasor suffers amplitude loss because the phasor is not 
completely “stopped.” (2) The signal will in general lie 
somewhere between Fourier resolution bins, so that the 
amplitude of the Fourier bin closest to the signal can be 
smaller than the true signal strength. The first effect can 
be ameliorated by discarding the Fourier resolution bins 
that lie toward the edges of the l/T Hz frequency span, at 
the cost of increasing the number of Fourier transforms re- 
quired to cover the entire search space of A v Hz. The frac- 
tion of Fourier resolution bins retained is referred to as \. 
The second effect can be reduced by “oversampling,” i.e., 
padding the data points with additional zeros, to decrease 
the frequency bin size of the Fourier transform. The de- 
gree of oversampling is described by a quantity, k, which is 
equal to N/(N + the number * of padded zeros). Note that 
decreasing k allows finer sampling of the Fourier transform 
at the cost of requiring increased computation (i.e., Fourier 
transforming a larger number of data points). 

A reasonable selection of \ and k, from the standpoint 
of minimizing acquisition time is \ = 1 and k = 1/2. For 
this choice of \ and « , it can be shown that the detection 
threshold parameter, £, is given by the following formula: 1 

1 S. A. Stephens, “An Analysis of FFT Tone Acquisition,” JPL In- 
teroffice Memorandum 335.1-92-14 (interned document), Jet Pro- 

pulsion Laboratory, Pasadena, California, May 14, 1992 


1/2 1/4 

0.995 = 2 J du J dw(SN R v ) 2 e-°^ aSNR ^ 

- 1/2 - 1/4 

CO 

x J dr r e- 0 5r2sNR ^I 0 (arSNR 2 v ) 
f 

where a = (sin (nu)/7ru) (sin ( ttw)/ww ), 7 0 is the 0th order 
modified Bessel function, and SNRy is th e voltage sign al- 
to-noise ratio over NT sec, \j2P t / N 0 NT . 

The probability that the noise amplitude in a given fre- 
quency bin exceeds £5 N Ry is exp (—0.5(£5'A^/^) 2 ). Thus, 
the probability, Pj a isei that the noise amplitude exceeds 
£SNR V in at least one of N frequency bins is 

Pfahe = 1 - (1 - exp ( -0.5(tSNR v ) 2 )) N 

The first equation allows numerical calculation of £ as 
a function of SNR V . Now if it is required that P/ a he = 
0.01, then the second equation allows calculation of N as 
a function of SNR V . But since SNRy = ^2P t /N 0 NT, N 
may then be calculated as a function of ( P t /No)T . Over 
the range 6 < SNR V < 12 } {P t / N 0 )T ^ 29N~° 92 . 

In the case of detection of a PN-modulated carrier, 
the calculation is analogous, but SNRy must be scaled 
by (1 — l/(2k c )) } where l/k c chips is the model sym- 
bol offset step size. This factor accounts for the SNRy 
loss during detection, due to the fact that the best align- 
ment of the model and signal sequences may be off by 
as much as l/2Ar c chips. If the search is conducted in 
half-chip increments ( k c = 2 ) and Pdetect = 0.99 while 
Pfahe = 0.01, then the relationship between ( P t /No)T 
and N is (P t /N 0 )T * 43iV" 0 - 91 
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